Magnonic spin-wave modes in CoFeB antidot lattices 
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In this manuscript time-resolved magneto-optical Kerr effect experiments on structured CoFeB films are pre- 
sented. The geometries considered are two dimensional square lattices of micrometer-sized antidots, fabricated 
by a focused ion beam. The spin-wave spectra of these magnonic crystals show a novel precessional mode, 
which can be related to a Bloch state at the zone boundary. Additionally, another magnetic mode of different na- 
ture appears, whose frequency displays no dependence on the externally applied magnetic field. These findings 
are interpreted as delocalized and localized modes, respectively. 

PACS numbers: 75.78.-n, 75.30.Ds, 75.50.Cc, 75.70.Ak, 75.40.Gb, 



Spin waves in confined magnetic structures have been in- 
vestigated for a few decades and the idea to consider peri- 
odically modulated structures can be dated back to 1976 \\\- 
Nowadays, magnetic metamaterials, i.e. materials with peri- 
odically varying magnetic properties in one, two or three di- 
mensions, are called magnonic crystals. The periodicity tai- 
lors the magnon band structure, eventually yielding band gaps 
and thus introducing the possibility to build (for example) 
spin-wave filters which have rather narrow transmission or re- 
jection bands . Ultimately, the control over spin waves is 
necessary for the long-term purpose of building magnetism- 
based logic devices that use spin waves rather than electrons 
for information processing. An important physical property 
of magnetic substructures is confinement caused by inhomo- 
geneous demagnetizing fields |5J- As a results, one finds lo- 
calized modes in magnonic crystals, which are restricted to 
regions of strong internal field inhomogeneity j6}-[8). These 
modes resemble atomic electronic states, which can extend 
throughout the whole crystal, if they get close enough to over- 
lap. Regarding magnonic crystals, this situation is accom- 
plished, if the inhomogeneity of the internal field connects the 
individual structural elements. Considering micrometer scales 
and an experimental approach in the time domain, the emer- 
gence and experimental proof of spin waves is suppressed 
when the magnetic damping is too high. Hence, for any ma- 
terial composing magnonic crystals, a small magnetic damp- 
ing is advantageous, as it is given for example in Permalloy 
(a = 0.008 [9|) or YIG (a = 0.00006 Q0)). In the present 
paper Co2oFe6oB2o has been investigated because of its par- 
ticularly small damping constant of a = 0.006 tl TTI . Also, 
CoFeB has the advantage of being amorphous, implying the 
absence of a magnetic anisotropy and grain boundaries and 
thus qualifying CoFeB as an ideal basic material for magnonic 
crystals. 

The samples were prepared by magnetron sputtering 
CoFeB films with a thickness d of 50nm onto a Si(100) sub- 
strate. A focussed ion beam technique (FIB) was used to 
create antidot structures with varying antidot radius r and 
lattice parameter a using a FEI Nova 600 NanoLab Dual- 
Beam. In all cases, the filling fraction was held constant at 



/ = nr 2 a~ 2 = 0.12, with the set of lattice parameters being 
3.5, 2.5 and 1.5 um. After the FIB processing a top layer of 
3 nm Ru was deposited by electron beam evaporation to pre- 
vent oxidation of the ferromagnetic layer. 

The experiments were carried out with an all-optical pump 
probe setup, which exploits in a time-resolved manner the 
magneto-optical Kerr effect (TRMOKE). High intensity pump 
pulses impinge on the sample surface with a fluence of 
57mJcm~ 2 and the successive temperature rise triggers the 
magnetic precession via an effective field pulse lfT2l[T3l . The 
excited high-k magnons decay within a few ps into spin waves 
with typical frequencies of a few GHz [14] and these mag- 
netic oscillations are detected with probe pulses, whose Kerr 
rotation upon reflection from the sample is proportional to the 
(dynamic) magnetization. With this technique, angular ampli- 
tudes of the excited spin waves of several degrees, depending 
on the pump fluence, can be reached. Note that the overall 
size of the structured regions was (150um) 2 . While this can- 
not be considered an infinitely extended crystal, it is still large 
enough for our experimental approach, as the pump (probe) 
pulses are focussed to a Gaussian width of 60pm (18|jm), re- 
spectively. The setup has previously been described in more 
detail elsewhere ||T5ll . 

Magnetization dynamics were measured covering a time in- 
terval of up to 1 ns after the arrival of the pump pulse. A quan- 
titative analysis was carried out by Fourier transformation. To 
minimize artefacts stemming from the (discrete) transforma- 
tion, a double exponential background originating from in- 
coherent magnons and phonons was subtracted from the raw 
data prior to calculation of the oscillations' FFT power. A 
series of TRMOKE spectra for systematically varied exter- 
nal fields can be processed this way and the resulting Fourier 
power spectra can be plotted in a color code in the frequency- 
field plane. In the present paper only readily analyzed data 
sets are presented, a more detailed example of the analysis 



is also given in Ref. [15]. From Fig. 1 a) it is apparent, that 
for small external fields the spin-wave amplitude strongly di- 
minishes, therefore normalization of each power spectrum to 
its respective maximum provides a better insight as shown in 



Fig. 1 b) The evolution of the observed magnetic modes with 
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the applied field can then be interpreted with respect to the 
origin of the different modes. 

The spin-wave theory for continuous films [16| provides 
the theoretical framework for the classification of the different 
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FIG. 1: (color online). Magnetization dynamics on a continuous 
CoFeB film of thickness d = 50 nm. (a) Fourier power spectrum 
as calculated from the time-resolved MOKE data 1151 . (b) After 
normalization the frequency branches can be attributed to the uni- 
form precession (white) as well as perpendicular standing spin waves 
(gray). Points are the peak positions, lines represent the fitted disper- 
sion {3} for n — and n = 1 , respectively. 
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FIG. 2: (color online). Not normalized (a) and normalized (b) power 
spectra of the magnetization dynamics on an array of antidots. The 
inset in (a) shows a SEM picture of an antidot structure and schematic 
of the field geometry. As compared to |Fig. 1| an additional magnetic 
mode is observed (orange), which can be identified as a Damon- 
Eshbach wave, owning the wavelength given by the periodicity of 
the antidots. 



modes. For an in-plane propagation direction which is per- 
pendicular to the applied field, spin waves are called Damon- 
Eshbach (DE) modes ifTTl . Neglecting the exchange contribu- 
tion, the dispersion relation valid for small k reads: 
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Within the same approximation the expression for an angle of 
45° between k and H reads: 
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If one allows for propagation in the direction normal to the 
film plane, one ends up with exchange dominated perpendic- 
ularly standing spin waves (PSSW): 
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For n = Eq. ([TJ describes the Kittel mode, which is charac- 
terized by a uniform precession of all spins. The abbreviations 
in Eqs. ([TJ, |2) and ^ are related to the external field H ext 
via (Oh = 7jUo#extCos30°, to the saturation magnetization Ms 
via (Om = YlkiMs, and to the exchange stiffness constant A via 

(Oa = Ym\ (f ) 2 ' wnere Y — g^B^ 1 is determined by the gy- 
romagnetic ratio. The trigonometric factor in (Oh accounts for 
the experimental necessity of a certain canting of the external 
field with respect to the sample's plane. 

The spectrum of the unstructured film shows two dif- 
ferent modes: the Kittel mode (« = 0) and the first order 
standing spin wave (n = 1). In combination with a SQUID 
measurement of Ms the following set of material parame- 
ters could be derived. The g-factor g = 2.04(1) as com- 
pared to 2.09 lfl8l is a little smaller than expected, whereas 
the saturation magnetization of jJ$Ms = 1.8(2) T and the 



low damping constant a = 0.006(1) match values found 
by other groups IfTTl [T9l . The exchange stiffness constant 
A = 2.0(1) x lO^'jirr 1 is smaller than the commonly used 
value 2.8(2) x 10" 11 Jm" 1 (20|, which belongs to the differ- 
ent composition Co72Fei8Bio- The Fourier power spectra in 
Fig.|2]show, that in the structured region a new mode appears, 
whose frequency is between the Kittel mode and the first order 
PSSW. The excitation mechanism does not select any propa- 
gation direction, wavelength or frequency. The population of 
the possible spin-wave modes is rather governed by the den- 
sity of states arising from the magnonic band structure. This 
density is high, wherever the dispersion is flat. According to 
this argument, a k- vector at the boundary of the first magnonic 
Brillouin zone can be expected, which means k = %cT 1 . As- 
suming a propagation direction perpendicular to the applied 
field, a fit of Eq. ([TJ to the newly appeared mode yields ap- 
proximately this value (see Fig. [2] and Fig.[3jd)). Higher or- 
der DE modes (k = nna^ 1 ) are less propable to detect due 
to damping. From the observation of the frequency alone one 
cannot judge, whether the propagation direction is 90° or 270° 
with respect to the external field, because the spin- wave theory 
ifTBIfTTIl provides both possibilities. The particular property of 
DE modes of being localized at one side of the film can be 
neglected in the present case, since In fact, it seems 

most probable, that spin waves travelling in both directions 
were excited. Then, the resulting superposition describes aid 
standing wave pattern, with vanishing group velocity. lETTl 

Rotating the sample in-plane to (p = 45° alters the spectrum, 
so that again only two modes can be identified (see Fig. At 
first sight the lower mode seems to be the Kittel mode as ob- 
served in the spectrum of the continuous film. A closer look 
reveals that this mode is shifted towards higher frequencies. 
The new mode for <p = 0° has a k-vector which points along 
the next nearest neighbour direction. For <p = 45° this crys- 
tallographic direction has a different orientation relative to the 
external field: now this mode should be described by Eq. ©. 
Fitting Eq. |2]i to the dominant peaks in Fig.[3ja) (yellow line) 
indeed yields again k = naT 1 . The measurements and analysis 
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FIG. 3: (color online), (a) to (c) Variation of antidot separation a for an external field orientation of <p = 45°. Solid yellow lines show fits of 
Eq. l|2j to field dependent modes (yellow dots), (d) Numerical computation of the internal field H ml = H ext +H^ em along the white section 
shown in the SEM picture. The external field IM)H ext = 140mT is oriented parallel to the section; the inset shows the wavelength as a function 
of the antidots distance a, stemming from the fits depicted in (a) to (c); orange dot from the fit of Eq. Qj, shown in Fig.[2jb). 



were carried out for all three different lattice constants. The 
experimental power spectra always show a mode, which is in 
frequency slightly above the Kittel mode and the fit confirmed 
in each case that k = Ka~ l (see Fig.^d)). For <p = 45° the ob- 
served quantized modes are much stronger than the quantized 
DE mode seen in the spectrum for <p = 0°. A simple rea- 
son could be, that the frequency is fourfold degenerate: spin 
waves with k = na^ 1 propagating in directions of 45°, 135°, 
225° and 315° relative to the external field have the same fre- 
quency and can be superposed to a 2d standing wave pattern. 
Although the discussed modes may possess stationary wave 
patterns, they originate from the delocalized modes described 
by the spin-wave theory of continuous films |[T6l . Owning 
their wavelength from the structural periodicity, these modes 
can be denoted as delocalized Bloch modes. 

It can be argued, that the detection of spin waves, which 
have a finite wavelength, with the given setup should at least 
be strongly suppressed by the fact, that the probe laser pulse 
averages over a certain area. Thus opposing phases contribute 
to the net signal. If they do so equally, the different contri- 
butions should indeed cancel each other. But since the probe 
pulses display a gaussian intensity profile in the lateral direc- 
tions, opposing phase contributions in different unit cells of 
the antidot lattice do not cancel totally. 

In addition to these findings, a second set of modes ap- 
pears with decreasing lattice constant, whose frequencies only 
weakly depend on the strength of the external field (see 
Fig. [3]). A simple ansatz to explain the field independence is 
to assume, that the observed oscillations belong to phononic 
excitations, such as surface acoustic waves. In the time- 
resolved reflectivity signal, indeed a weak oscillation is ob- 
served. The associated frequency of roughly 0.5 GHz does not 
compare to the frequencies which appear in the time-resolved 
MOKE spectra. Although the discussed modes do not sig- 
nificantly depend in frequency on the strength of the applied 
field, changing the relative orientation of i/ ext alters this fre- 
quency, or introduces a second or third mode with a fixed but 



different frequency. Also, resonance effects can be seen in 
not normalized spectra, when the frequencies of the different 
sets of modes intersect. This clearly rules out a phononic ori- 
gin. The following consideration relates the appearance of the 
new mode to the internal field and suggests that the mode is 
a localized spin wave, existing close to each antidot's borders 
and oscillating in phase throughout the whole crystal. The 
geometric filling fraction / = 0.12 remains the same for the 
different lattice constants, whereas the fraction of the regions 
with a not saturated, inhomogeneous internal field grows with 
a shrinking lattice constant. A numerical computation of the 
internal field, using the software package Nmag [22], illus- 
trates this behavior. The internal field is given in Fig. [3jd). 
We expect that localized spin-wave modes should therefore 
become more and more dominant. For a = 3.5pm the new 
modes can almost not be distinguished from noise and the 
modes known from the continuous film and the delocalized 
mode dominate the spectrum. For a — 1 .5 um, this situation 
is completely reversed, as can be seen in Fig. |3jc). Now re- 
gions with a saturated internal field H mt — H ext only become 
relevant enough in size to significantly allow for delocalized 
spin waves, when the external field exceeds a certain threshold 
value, which increases with a decreasing lattice constant. Il23l 
This shows, that apart from the absolute size of the structures, 
the external field can be used to tune the spin-wave spectrum 
of the same magnonic crystal between a metallic and a more 
atomic behavior, with delocalized and localized spin waves 
dominating, respectively. 

In summary an all-optical pump-probe technique was used 
to investigate spin-wave spectra of structured CoFeB films. 
The spectra of this kind of 2d magnonic crystal show modes 
with a k- vector k = naT x induced by the period a of the 
magnonic crystal. Reducing the lattice constant from 3.5 to 
1.5um, new modes apear which dominate for small a. Our 
explanation suggests, that these modes are localized in re- 
gions of strongly inhomogeneous internal field. Research was 
funded by DFG priority program 1 133. 



4 



[1] C. Sykes, J. Adam, and J. Collins, Applied Physics Letters 29, 
388 (1976). 

[2] M. Krawczyk and H. Puszkarski, Physical Review B 77, 
054437 (2008). 

[3] A. V. Chumak, A. A. Serga, B. Hillebrands, and M. P. Kostylev, 
Applied Physi cs Letters 93, 022508 (2008), URL |http://| 
| link . aip . org/link/?APL/93/022508/T1 

[4] K.-S. Lee, D.-S. Han, and S.-K. Kim, Phys. Rev. Lett. 102, 
127202 (2009). 

[5] R. Joseph and E. Schlomann, Journal of Applied Physics 36, 
1579 (1965). 

[6] S. Neusser, B. Botters, M. Becherer, D. Schmitt-Landsiedel, 
and D. Grundler, Applied Physics Letters 93, 122501 ( 2008), 
URL [http : //link . aip . org/link/?APL/93/12250l7l1 
[7] J. Jorzick, S. O. Demokritov, B. Hillebrands, M. Bailleul, 
C. Fermon, K. Y. Guslienko, A. N. Slavin, D. V. Berkov, and 
N. L. Gorn, Phys. Rev. Lett. 88, 047204 (2002). 
[8] M. Pechan, C. Yu, R. Compton, J. Park, and P. Cowell, Journal 

of Applied Physics 97, 10J903 (2005). 
[9] J. Walowski, G. Miiller, M. Djordjevic, M. Munzenberg, 
M. Kliiui, C. A. F. Vaz, and J. A. C. Bland, Phys. Rev. Lett. 
101, 237401 (2008). 

[10] J. Goulon, A. Rogalev, F. Wilhelm, N. Jaouen, C. Goulon- 
Ginet, G. Goujon, J. B. Youssef, and M. Indenbom, Journal of 
Electron Spectroscopy and Related Phenomena 156-158, 38 
(2007), ISSN 0368-2048, electronic Spectroscopy and Struc- 
ture: ICESS-10, URL http : //www. sciencedirect . 

| com/science/article/B6TGC-4MG5WB0-2/2/ | 
259da6be81el612c7e421a3f221dedf9 

[11] M. Oogane, T. Wakitani, S. Yakata, R. Yilgin, Y. Ando, 
A. Sakuma, and T. Miyazaki, Japanese Journal of Applied 
Physics 45, 3889 (2006). 

[12] M. van Kampen, C. Jozsa, J. T. Kohlhepp, P. LeClair, L. Lagae, 
W. J. M. de Jonge, and B. Koopmans, Physical Review Letters 
88, 227201 ( 2002), URL |http : //link. aps . org/abstract/| 
PRL/v88/e22720H 

[13] Y. Liu, L. R. Shelford, V. V. Kruglyak, R. J. Hicken, 



Y. Sakuraba, M. Oogane, Y. Ando, and T. Miyazaki (2007), vol. 
101, p. 09C106, URL |http : //link . aip . org/link/? JAP/| 
101/09C106/1I 

[14] M. Djordjevic and M. Munzenberg, Physical Review B 75, 
012404 (2007), U RL |http : //link . aps . org/abstract7| 
PRB/v75/e012404l 

[15] J. Walowski, M. D. Kaufmann, B. Lenk, C. Hamann, J. Mc- 
Cord, and M. Munzenberg, Journal of Physics D: Applied 
Physics 41, 1 64016 (2008), URL http : //stacks . iop . org/| 
0022-3727/41/i=16/a=16 4016 

[16] B. Kalinikos and A. Slavin, Journal of Physics C: Solid State 
Physics 19, 7013 (1986). 

[17] R. Damon and J. Eshbach, Journal of Physics and Chemistry of 
Solids 19, 308 (1961). 

[18] P. Paluskar, R. Lavrijsen, M. Sicot, J. Kolhepp, H. Swagten, and 
B. Koopmans, Physical Review Letters 102, 016602 (2009). 

[19] R. O'Handley, R. Hasegawa, R. Ray, and C.-P. Chou, Applied 
Physics Letters 29, 330 (1976). 

[20] C. Bilzer, T. Devolder, J. Kim, G. Counil, and C. Chappert, 
Journal of Applied Physics 100, 053903 (2006). 

[21] In comparison with the unstructured film the Kittel mode seems 
much stronger in power inside the structure. Note that a quan- 
tized Backward- Volume mode (k||H) with k = 7la~ l has a fre- 
quency, which is so close to the Kittelmode, that the frequency 
resolution in the experiment (0.5 GHz) is too large to distin- 
guish the modes. This consideration indicates, that quantized 
Backward- Volume modes may also contribute to the oscilla- 
tions denoted simply as uniform precession in Fig.[2]. 

[22] T. Fischbacher, M. Franchin, G. Bordignon, and H. Fangohr, 
IEEE Transactions on Magnetics 43, 2896 (2007). 

[23] This can also be seen in numerical simulations of the internal 
field. To avoid overcrowding, a field variation is not included 
in Fig.[3fd). Considering the fraction of the inhomogeneous in- 
ternal field, the effect of reducing the lattice constant is qual- 
itatively the same as for decreasing the external field: in both 
cases the fractional area of non-saturated H mt grows.. 



